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ABSTRACT 



Aims. We aim to investigate the relation between the long-term flux density and the position angle (PA) evolution of inner-jet in 
blazars. 

Methods. We have carried out the elliptic Gaussian model-fit to the 'core' of 50 blazars from 15 GHz VLBA data, and analyzed the 
variability properties of three blazars from the model-fit results. 

Results. Diverse correlations between the long-term peak flux density and the PA evolution of the major axis of the 'core' have been 
found in ~ 20% of the 50 sources. Of them, three typical blazars have been analyzed, which also show quasi-periodic flux variations 
of a few years (T). The correlation between the peak flux density and the PA of inner-jet is positive for S5 0716+714, and negative 
for S4 1807+698. The two sources cannot be explained with the ballistic jet models, the non-ballistic models have been analyzed to 
explain the two sub-luminal blazars. A correlation between the peak flux density and the PA (with a T/4 time lag) of inner-jet is found 
in [HB89] 1823+568, this correlation can be explained with a ballistic precession jet model. All the explanations are based mainly on 
the geometric beaming effect; physical flux density variations from the jet base would be considered for more complicated situations 
in future, which could account for the no or less significance of the correlation between the peak flux density and the PA of inner-jet 
in the majority blazars of our sample. 
Conclusions. 

Key words. BL Lacertae objects: individual: S5 0716+714, 3C 371, 4C +56.27 - radio continuum: galaxies - galaxies: jets 
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1. Introduction 

Blazars are thought of the extreme of active galactic nuclei 
(AGNs), which consist of flat-spectrum radio quasars and BL 
Lac objects. They are variable in almost whole electromagnetic 
spectrum from radio to gamma-ray. The extreme properties of 
blazars are mainly attributed to their jets which are pointing 
closely to our line of sight, according to the unified model of 
AGN (Urry & Padovani 1995). With very long baseline interfer- 
ometry (VLBI) - the highest angular resolution to date in astron- 
omy - blazars can be resolved into a bright core and often one- 
sided jet components in parsec scale or sub-parsec scale (Lister 
et al. 2009a). In the VLBI scales, the blazar jets frequently 
show fast outward motions, with an apparent advancing speed 
of greater than the speed of light - namely super-luminal mo- 
tion (Rees 1966). This is due to the so-called relativistic Doppler 
boosting - relativistic beaming effect (e.g. Cohen et al. 2007; 
Hovatta et al. 2009). Such high flux variability and VLBI struc- 
ture changes together provide us an effective means to study the 
inner-jet property and probably the central engine's property of 
the blazars in further. 

Recently, Britzen et al. (2009) found there seems to be evi- 
dence for an apparent stationarity of jet components (with regard 
to their core separation) with time, in blazar S5 0716+714, while 
the inner-jet components exhibit strong changes in their position 
angles. This indicates that the jet components move non-radially 
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(or non-ballistic) with regard to their position angle. These au- 
thors attributed such sub-luminal motions to a geometric ori- 
gin, although the geometric models were not well defined yet. 
Furthermore, they found that the long-term flux density and the 
position angle of the jet in 0716+714 show a significant positive 
correlation. The similar phenomena of the apparent stationar- 
ity of jet components with time have also been reported in BL 
Lac objects S5 1803+784 and PKS 0735+178, by Britzen et al. 
(2010a, 2010b). From CJF sample (Caltech-Jodrell Bank flat- 
spectrum sample of radio loud active galaxies), Karouzos et al. 
(2012) investigated the morphologies and pc-scale jet kinemat- 
ics of more than 200 sources by using the so-called jet ridge-line 
method. They found that about a half of the sample show jet 
widths > 10°, with BL Lac jet ridge lines showing significantly 
larger apparent widths than both quasars and radio galaxies, in- 
dicating that the BL Lacs have more significant position angle 
changes of jets in their evolution. 

From the study of AGN jet kinematics both in the 2 cm sur- 
vey (Kellermann et al. 2004) and in the radio reference frame 
image database (RRFID), about one-third of well-measured 
component trajectories are non-radial (Piner et al. 2007). In 
the MOJAVE (monitoring of jets in active galactic nuclei with 
VLBA experiments) sample, 166 out of the 526 robust compo- 
nents (32%) show non-radial motions, indicating that non-radial 
motion is a common feature of the jet flow in blazars (Lister et al. 
2009a, 2009b). Three of these non-radial components might be 
considered inward at the 3cr level (Lister et al. 2009b), this could 
be due to a projection effect of curved jets. With the acceleration 
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measurements of the components of MOJAVE sources, Homan 
et al. (2009) found the perpendicular accelerations are closely 
linked with non-radial motions, and about half of the compo- 
nents show 'non-radial' motions; but the parallel accelerations 
are generally larger than the perpendicular accelerations with re- 
spect to the component velocities. 

To further investigate the relation, which was first discovered 
in 0716+714 by Britzen et al. (2009), between the long-term 
flux density and the position angle (PA) evolution of inner-jet 
in blazars, we model-fitted the 'core' components of 50 core- 
dominated balzars from the MOJAVE data released up to 201 1 
(Lister et al. 2009a). We aimed at finding more sources which 
showing correlations between the long-term flux density and the 
PA evolution of inner-jet. We also analyzed the inner-jet models 
to explain the typical correlations, if any, in the model-fit results. 

2. Model-fit result and analysis 

2.1. Model-fit to source 'core' 

We selected 50 core-dominated blazars showing one-sided VLBI 
core-jet and they have been observed with the very long base- 
line array (VLB A) at 15 GHz in more than 10 years, from 
the MOJAVE database (Lister et al. 2OO9£0). We model-fitted 
the 'core' component of each source, with an elliptic Gaussian 
model in the AIPS task 'JMFIT', to get the peak flux density 
and the PA of the major axis of the Gaussian component, etc. 
We consider that the 'core' in the 15 GHz VLBA images of these 
core-dominated blazars is the 'inner-jet' rather than the true core. 
Therefore, the inner-jet can be modeled with an elliptic Gaussian 
component which its major axis is along with the inner-jet ori- 
entation or the inner-jet ridge-line on average, thus reflecting the 
inner-jet PA. 

We are interested in the long-term flux density and the PA 
evolution of the 'core' of blazars, 15 GHz MOJAVE datasets 
are the best because these observations have been spanned over 
10 years. With the calibrated visibility of the MOJAVE data we 
could do the model-fit to whole source structure, however, for the 
50 sources with more than 10 epochs of each source we couldn't 
complete the multi-components model-fit to all the source struc- 
tures, as usually did with the modelfit program in the Difmap 
package. We notice that it is not needed currently to model-fit 
whole source structure, for that our purpose is studying on aver- 
age the innermost jet behavior of the sources. An advantage of 
our approach is that we could get ride of an ambiguity in iden- 
tifying VLBI components from epoch to epoch which have long 
been discussed as a potential problem for multi-epoch VLBI ob- 
servations (see discussions in e.g., Vermeulen et al. 2003; Piner 
et al. 2007; Britzen et al. 2009). Our approach is similar to that 
the model-fitting to the 'core' of 250 flat spectrum radio sources 
by Kovalev et al. (2005) with the visibility function. 

The solution of the JMFIT in AIPS gives us the peak flux 
density per beam, the X and Y position of a Gaussian compo- 
nent, the de-convoluted major/minor axis and the PA of the ma- 
jor axis of the Gaussian component. The internal error in the 
JMFIT is very small. We have done twice the JMFIT for each 
dataset in order to get a difference as the error of the solution. 
The final error of each fitting parameter we adopted, is the larger 
one between the internal fitting error and the difference of the 
twice fittings. 

We have obtained the model-fit results of more than 500 
epochs, we here report only for three sources that have shown 

1 http://www.physics.purdue.edu/MOJAVE/ 



typical correlations between the long-term flux density and the 
PA of the major axis of 'core'. We present the model-fit results 
for the three sources: namely S5 0716+714, S4 1807+698, and 
[HB89] 1823+568 in TableQ] The columns give; (1) the source 
name; (2), (3), (4), (5), (6), (7) the minimum /median /maximum 
value of peak flux density per beam of elliptic Gaussian compo- 
nent, major axis of the component, minor axis of the component, 
position angle of the major axis, component position shift from 
the reference point [0,0], position angle of the component posi- 
tion shift, respectively; (8) linear Pearson correlation coefficient 
(and significance) between the peak flux density and the position 
angle of the major axis; (9) linear Pearson correlation coefficient 
(and significance) between the position angle of major axis and 
the position angle of component shift; (10) major period of the 
peak flux variations derived from the power spectrum density 
analysis. 

It should be noted that there are no correlations between the 
parameters of the de-convoluted elliptic Gaussian model from 
the model-fit and the parameters of the restoring beam of im- 
age. For instance, there is no any correlation between the PA of 
the major axis of the 15 GHz VLBA 'core' and the PA of the 
restoring beam of the 15 GHz VLBA image for 0716+714, as 
shown in Fig.Q] indicating that the model-fit result is related to 
the source-intrinsic property rather than related to the restoring 
beam of image. 

2.2. Correlation and PSD analysis 

The linear Pearson correlation is analyzed in the model-fit pa- 
rameters. A diversity of correlations between the long-term peak 
flux density and the PA evolution of the major axis of the 'core' 
has been found in ~ 20% of 50 sources, including positive cor- 
relation, anti-correlation, and correlation with some time delay. 
There shows the positive, the negative and the time-delayed cor- 
relation between the peak flux density and the PA of the ma- 
jor axis for blazars 0716+714, 1807+698 and 1823+568 respec- 
tively in TableQ](also see, Fig. [2] Fig. [3] and Fig. More dis- 
cussions on the correlations are given in the section 3. However, 
nearly 80% of the 50 sources have shown no or low significance 
of correlation between the long-term peak flux and the PA evo- 
lution of major axis. 

We also analyzed the possible periodicity of the long-term 
peak flux variations, with the power spectrum density (PSD) 
method for the three sources. All the three sources show a 
peak in the PSD plot (Fig. [5]), the resulted major periodicity 
for 0716+714, 1807+698 and 1823+568, is respectively listed 
in Table [TJ 

2.3. The three typical sources 
2.3.1. S5 0716+714 

S5 0716+714 is a well known BL Lac object, with a redshift of 
0.31 (Nilsson et al. 2008). The source is variable from radio to 
gamma-ray (Wagner et al. 1996; Raiteri et al. 2003; Liu et al. 
2012; Zhang 2010; Abdo et al. 2010), and it is also a TeV source 
(Tavecchio et al. 2010). The VLA map at 20 cm shows kpc-scale 
jets embedded in an oval cocoon (Antonucci et al. 1986). VLBI 
images exhibit an extremely core-dominated jet pointing to the 
north and misaligned with the VLA jet by ~ 90°, e.g., at 5 and 
8.4/22 GHz (Gabuzda et al. 1998). The 15 GHz MOJAVE im- 
ages show a bright compact core and a weak jet within 5 mas, 
with a jet advancing speed of 10 c (Lister et al. 2009b). However, 
the cross-epoch identification of components is not unique, from 
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Table 1. The model-fit result of inner-jet in three blazars from the MOJAVE 15 GHz data, the result of correlation analysis, and the 
major period of peak flux variations derived from the PSD. 
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Source 


Peak flux 


Major axis 


Minor axis 


PAofM.A. 


Dist. from ref. 


PA of comp. 


Corr.2&5 


Corr.5&7 


T 




[Jy/beam] 


[mas] 


[mas] 


[deg] 


[mas] 


[deg] 






[yi 


0716+714 


0.2258/1.4631/3.3709 


0.11/0.21/0.37 


0.02/0.05/0.09 


6.6/20.6/43.0 


0/0.01/0.09 


-7.3/20.5/50.4 


0.44(4.6E-4) 


0.62(2. 1E-7) 


5.8 ± 


1807+698 


0.4826/0.6987/0.9561 


0.50/0.68/0.91 


0.09/0.11/0.23 


-110/-105/-99 


0.04/0.07/0.14 


-112/-104/-97 


-0.45(0.02) 


0.21(0.27) 


8.3 ± 


1823+568 


0.5676/1.1359/2.0515 


0.31/0.40/0.70 


0.02/0.06/0.13 


-165/-161/-151 


0.02/0.03/0.08 


-147/-146/-142 


0.56(1. 6E-5)° 


0.29(0.03)° 


7.0 ± 



Notes. "Epoch of the PA of the major axis in 1823+568 data has been shifted forward by T/4 in the correlation analysis. 



some different data Britzen et al. (2009) identified the jet com- 
ponents which are almost stationary respect to the core. 

The correlation between the peak flux density and the PA of 
the major axis of the 'core' in 0716+714 is significant (Table[TJi. 
Though during the first part of the double-peaks in 2009-2010, 
the PA is not significantly changed but it does follow the trend 
of flux density variations (see Fig. |2j. From long timescale radio 
data, a 5.7+0.5 yr period is reported (Raiteri et al. 2003; Fan et al. 
2007). From our model-fit data at 15 GHz, the peak flux density 
does show a peak in 2003-2004 and double peaks in 2009-2010 
(Fig. with an overall periodic flux variations of 5.8+0.4 yr 
in the PSD plot (Fig. [5] only data after 2001.0 were used in the 
PSD, because the VLBA data before 2001.0 are poor). There is 
also a possible period of 1.2+0.2 yr in Fig. [5] 

2.3.2. S4 1807+698 

The source 1807+698 (3C 371, z=0.051) is related to BL Lac 
objects because of its optical variability and polarization in a 
large elliptical galaxy whose emission is dominated by the flux 
from the nucleus (Fan et al. 1999; Katajainen et al. 2000), and it 
is detected in gamma-ray by Fermi LAT. VLA image at 1.36 
GHz shows a kpc scale halo and jets elongated in east-west 
(Cassaro et al. 1999). VLBI maps reveal a one-sided core-jet 
at PA ~ -100°, e.g., at 5 GHz (Gomez & Marscher 2000). The 
jets in the optical with HST and in the X-rays with Chandra have 
been detected (Pesce et al. 2001), which are roughly aligned with 
the parsec-scale radio jet. The 15 GHz MOJAVE result shows a 
sub-luminal jet motion of 0.1 c (Lister et al. 2009b). A 43 GHz 
VLBA image exhibits a jet flow assembling a wiggling pattern 
(Lister 2001). 

An anti-correlation between the peak flux density and the PA 
of the major axis of inner-jet, with the coefficient of -0.45 at a 
confidence level of 98%, is found, and a flux variability period 
of 8.3+0.6 yr is estimated with the PSD method (see, Fig.[3]and 
Fig.©. 



2.3.3. [HB89] 1823+568 

The source 1823+568 (4C +56.27, z=0.664), is also a BL Lac 
object, and it is detected in gamma-ray by Fermi LAT. It shows a 
halo and one-sided kpc-scale jet to east (Kollgaard et al. 1992). 
VLA observation at 5 GHz by O'Dea, Barvainis & Challis 
(1986) shows a curved jet extending southward from the core 
and turning to the east. The object is highly variable and has a 
strong linear polarization both in radio and optical bands (Lister 
2001). VLBI maps reveal a compact core and jets extending to 
the south, e.g., at 5 GHz (Fomalont et al. 2000). The 15 GHz 
MOJAVE result shows a super-luminal jet motion of 20.8 c 
(Lister et al. 2009b). 




-60 -40 -20 20 40 
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Fig. 1. PA of the major axis of the 15 GHz VLBA 'core' (inner- 
jet) versus PA of the restoring beam of the 15 GHz VLBA image 
for 0716+714. 
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Fig. 2. Peak flux density and position angle of the major axis of 
inner-jet in 0716+714 from the MOJAVE 15 GHz data versus 
epoch. Vertical lines show the estimated intervals from peak to 
peak of flux, and compared with the result from PSD analysis. 



Our model-fit result from the MOJAVE 15 GHz 'core', sug- 
gests a significant correlation between the peak flux density and 
the T/4 time-shift-forward PA of the major axis of inner-jet in 
TableQ] The variability period of T=7.0+0.4 yr is estimated with 
the PSD (see, Fig. Hand Fig. 0. 



3. Inner-jet models and discussion 

We summarize our findings from the three BL Lac objects. 
There exhibit diverse correlations between the peak flux density 
and the PA of the major axis of 'core' component in the three 
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Fig. 3. Peak flux density and position angle of the major axis of 
inner-jet in 1807+698 from the MOJAVE 15 GHz data versus 
epoch. Vertical lines show the estimated intervals from peak to 
peak of flux, and compared with the result from PSD analysis. 
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Fig. 4. Peak flux density and position angle (shifted forward by 
1.7 yr) of the major axis of inner-jet in 1823+568 from the 
MOJAVE 15 GHz data versus epoch. Vertical lines show the es- 
timated intervals from peak to peak of flux, and compared with 
the result from PSD analysis. 



sources. They also show quasi-periodic variations of the peak 
flux density, around a few years. 

To explain the correlations and periodicity mentioned above, 
we have to find the models related to the inner-jet process in 
the radio loud AGNs. In general, from the literature there are 
several parameters/factors related to jet models, which should 
be clarified and clearly defined, e.g., 'ballistic' (we define as a 
single jet component always moving along a straight line once 
launched from the jet nozzle), 'non-ballistic' (we define as a sin- 
gle jet component moving in a curved way), 'precession' (this 
has been confused in the literature, we define it as a precession 
'jet' in terms of a jet flow (by no means a single jet component), 
in which every single jet component is ballistically launched 
from a precession jet-nozzle), 'non-precession' (we define as a 
jet nozzle has no precession). There are four meaningful hybrids 
of the parameters mentioned above, i.e., ballistic jets with a non- 
precession nozzle (B+nP), ballistic jets with a precession nozzle 
(B+P), non-ballistic jets with a non-precession nozzle (nB+nP), 
non-ballistic jets with a precession nozzle (nB+P). The quasi- 
periodic changing of the PA of inner-jet in the three sources, 
could be explained with the hybrid models except the first one 
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Fig. 5. Power spectrum density of the peak flux density of 
0716+714 (data with epochs > 2001.0 were used, upper panel), 
1807+698 (middle panel) and 1823+568 (lower panel), respec- 
tively. 




Fig. 6. A schematic diagram of the ballistic jets with a preces- 
sion nozzle (i.e. the 'B+P' model), with a counter-clockwise pre- 
cession 'B— >A— >D— »C' for 1823+568, a jet at the phase 'A' is 
pointing the closest to our line of sight. 



0716+714 



1807+698 





Fig. 7. A schematic diagram of a non-ballistic jet with a counter- 
clockwise movement A->B->C->D' for 0716+714 (left), and a 
non-ballistic jet with a clockwise movement A— >B— >C— »D' for 
1807+698 (right). In both cases, a jet at phase 'D' is pointing the 
closest to our line of sight. 



4 



X. Liu et al.: Correlations between the peak flux density and the position angle of inner-jet in three blazars 



(B+nP). A precession of the jet nozzle can be due to e.g., dou- 
ble super-massive black holes or other reasons (see Britzen et al. 
2010a and references therein). For the non-ballistic jet, models 
have been proposed e.g., Gong et al. (201 1), also see Lister et al. 
(2009b) and Britzen et al. (2009, 2010a, 2010b) for discussions. 
In these models, a non-ballistic motion of jet is largely related 
to the jet interaction with ambient matter. In the following, we 
will use the hybrid inner-jet models we classified, for the three 
sources. 

Firstly, we try to apply these models to 1823+568, after anal- 
ysis below, the B+P model (i.e. ballistic jets with a precessionjet 
nozzle) is applicable to this source, as shown in the schematic 
diagram of Fig. [6] We are nearly face-on the jet with the esti- 
mated jet viewing angle of ~8.4°(Savolainen et al. 2009), and the 
inner-jet PA is oscillating between -151°and -165°(see TableQ] 
and Fig.|4|. In Fig. [6] our line of sight is close to the phase A' 
of a precessionjet flow. At the phase A', the relativistic beam- 
ing effect is the strongest, so the flux density peaks; on the con- 
trast, the phase 'C of precessionjet flow is far from our line of 
sight, so the flux density is the lowest; at the phase 'B' or 'D' 
the flux density has a middle value due to the beaming effect. 
For the periodic peak flux variations (T~7.0 yr in Table Q]), and 
the correlation between the peak flux density and the PA (with a 
T/4 time-lag) of inner-jet, the precession of the jet flow must be 
counter-clockwise, i.e. through the phases 'B— >A— >D— >C— >B' 
in Fig. [6] Such that the PA peaks (a minus value) at the phase 
'B' while the flux density is at a middle level; and the PA has a 
middle value at the phase A' while the flux density peaks; and 
the PA has the smallest (minus) value at the phase 'D' while the 
flux density is at a middle level; and the PA has a middle value 
at the phase 'C while the flux density is the lowest (also see, 
Fig- HI- So the correlation between the peak flux density and the 
PA (with the T/4 time lag) of inner-jet can be well explained. 
This 'B+P' model seems to support the super-luminal motions 
of 20.8 c (Lister et al. 2009b) or 9.4 c (Savolainen et al. 2009) in 
this source, because the ballistic motions have generally faster 
radial speed than that of the non-ballistic motions. The origin of 
the periodic precession of ~7.0 yr in the jet nozzle still needs to 
be investigated in future. 

For the source either S5 0716+714 or 1807+698, the cor- 
relations between the peak flux density and the PA of inner-jet 
cannot be explained with a ballistic jet model. For 0716+714, 
the inner-jet PA is oscillating between 6.6°and 43°(see Table [1] 
and Fig. [2]). We are nearly face-on the jet with the estimated 
jet viewing angle of ~5.3°(Savolainen et al. 2009). The signif- 
icant positive correlation (Table [TJ indicates that the flux den- 
sity peaks while the PA peaks (i.e. far from North) simultane- 
ously (see Fig. |2|, so it is impossible that the jet at the phase 
A' in the schematic diagram (Fig. [7] left) has a middle-level of 
flux density according to the beaming effect within a ballistic 
model. Non-ballistic jet models have to be invoked, in which 
a single jet component is moving in a curved way. In the non- 
ballistic scenario, regardless the jet nozzle is in precession or 
not, the jet must have a non-radial motion. For 0716+714, in the 
schematic diagram Fig. [7] (left), for a non-ballistic jet moving 
along the counter-clockwise helical trajectory A— >B— >C— >D', 
i.e. the curved way through which the jet at the phase 'D' is 
pointing the closest to our line of sight, such that the flux density 
peaks while the PA peaks too. The peak flux density has a mid- 
dle value while the PA has a middle value too at the phase A', 
and so on, the positive correlation between the peak flux density 
and the PA at other phases 'B', 'C can also be explained with 
the counter-clockwise non-ballistic helical jet model (see Fig. 2 
and Fig. [7] left). The non-ballistic jet model will lead to a less 



outward motion speed than that of a ballistic jet model, assum- 
ing the originally-launched jet speed is the same. This is con- 
sisting with the sub-luminal motion or the 'stationary scenario' 
of 0716+714, as proposed by Britzen et al. (2009). In general, a 
non-ballistic jet is likely to be formed, through the shocked jet 
plasma working on the surrounding medium (Gong et al. 201 1), 
and a helical like trajectory of the jet ridge-line would be ex- 
pected. 

For 1807+698, the inner-jet PA is oscillating between 
-99°and -1 10°(see TableQ]and Fig. [3). We are not well face-on 
the jet for the estimated jet viewing angle of ~45.3°(Savolainen 
et al. 2009). A clockwise non-ballistic jet motion can explain 
the anti-correlation (i.e. a T/2 time lag) between the peak flux 
density and the PA of inner-jet, i.e. a jet is moving through 
A->B->C->D' in Fig.|7]right, the jet through at the 'D' is point- 
ing the closest to our line of sight. Such that the flux density 
peaks due to beaming effect at the 'D' while the PA has a small- 
est (minus) value, and at the phases A' and 'C', both the peak 
flux density and the PA have middle values, and at the phase 'B' 
the flux density is the lowest while the PA has the largest (minus) 
value. The source jet is sub-luminal (0. 1 c, Lister et al. 2009b), 
which is the non-ballistic jet model expected. The jet viewing 
angle of 1807+698 is relatively large (45.3°), so that the beam- 
ing effect is reduced. This may have led to the relatively low 
significance of the correlation between the peak flux density and 
the PA of inner-jet of this source (see Table Q]), and the lower 
significance of the estimated period of the flux variations. 

Among the three sources, the periodicity of flux variabil- 
ity and the correlation between the peak flux density and the 
T/4-laged PA of inner-jet in 1823+568 is a purely geometric ef- 
fect according to the ballistic precessionjet model 'B+P'. The 
quasi-periodicity of flux variability and the correlations between 
the peak flux density and the PA of inner-jet in 0716+714 and 
1807+698, as the non-ballistic models considered above, are 
also mostly caused by the geometric beaming effect. For the non- 
ballistic models in particular, however, even if the precession of 
the jet nozzle is invoked, the repeating periodic flux variability 
frequently found in blazars still needs to be fully explained. A 
kind of physical (non-geometric) outburst from the jet nozzle, 
or a physically periodic flux variations from the jet base would 
be considered for 0716+714 and 1807+698 in the future, which 
may be related to perturbations induced by the changes of the 
accretion mode or by a periodic intervening of a companying 
object. Multi-bands study of the three blazars on e.g., the SED, 
is also needed, like that recently being investigated for balzars 
statistically (see, Meyer et al. 201 1; Lister et al. 201 1). 

4. Summary 

We have carried out the elliptic Gaussian model-fit to the 'core' 
of 50 blazars. Diverse correlations between the long-term peak 
flux density and the PA evolution of the major axis of the 'core' 
have been found in ~ 20% of the 50 sources. Of them, three 
typical blazars have been analyzed, which also show quasi- 
periodic flux variations of a few years (T). The correlation be- 
tween the peak flux density and the PA of inner-jet is positive for 
0716+714, and negative for 1807+698. The two sources cannot 
be explained with the ballistic jet models, the non-ballistic mod- 
els have been analyzed to explain the two sub-luminal blazars. A 
correlation between the peak flux density and the PA (with a T/4 
time lag) of inner-jet is found in 1823+568, this correlation can 
be explained by a ballistic precessionjet model. All the explana- 
tions are based mainly on the geometric beaming effect; physical 
flux density variations from the jet base would be considered for 
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more complicated situations in future, which could account for 
the no or less significance of the correlation between the peak 
flux density and the PA of inner-jet in the majority blazars of our 
sample. 
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